Solitary waves in elongated clouds of strongly-interacting bosons 
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We examine the propagation of solitary waves in elongated clouds of trapped bosonic atoms as 
the confinement, the strength of the interatomic interaction, and the atom density are varied. We 
identify three different physical regimes and develop a general formalism that allows us to interpolate 
between them. Finally we pay special attention to the transition to the Tonks-Girardeau limit of 
strongly-interacting bosons. 
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By appropriate manipulation of the trapping poten- 
tial, confined ultracold atoms can achieve conditions of 
reduced dimensionality. Taking advantage of this prop- 
erty, two recent experiments described in Refs. 0,0 have 
managed to create conditions such that bosonic atoms 
approached the so-called Tonks-Girardeau phase |3| pre- 
dicted for strongly-interacting bosons in one dimension. 
In the experiment of Ref. [2j , bosonic atoms were confined 
in elongated traps, and the gas was observed to approach 
the Tonks-Girardeau limit as the transverse confinement 
was increased. Evidence for this transition was obtained 
from measurements of the energy and of the axial size 
of the gas. In the experiment of Ref. Q the momentum 
distribution provided evidence for this transition. 

Motivated by these experiments, we will consider here 
how solitary waves emerge in strongly-interacting, quasi- 
one-dimensional atomic gases. We first present a gen- 
eral description of the problem, identifying three phys- 
ically distinct regimes |^y. We then develop a formal- 
ism that allows us to calculate the density profile, en- 
ergy, and momentum associated with grey/dark solitary 

waves HHHIEIIHIIIIII- This formalism can be 
applied to any quasi-one dimensional system given only 
knowledge of the energy per unit length as function of 
the density per unit length. We find interesting changes 
in the properties of solitary waves as a function of density 
which suggest that the study of solitary waves in these 
elongated clouds could provide confirmation of the tran- 
sition to the Tonks-Girardeau limit 0, H3| • Our results 
also suggest eventual technological applications, e.g., the 
transmission of signals in atomic waveguides. 



the form 



ihd t V = (-h 2 \/ 2 /2M + U \^\ 2 + V)^. 



(1) 



Here, we will neglect the effects of trapping along the z- 
axis of weak confinement and will for simplicity consider 
a cylindrical trap, V = Mlo\(x 2 + y 2 )/2, where M is 
the atom mass, and to± is the frequency of the trapping 
potential in the transverse direction. We approximate 
a short-ranged atom-atom interaction as VS n t(r — r') = 
UqS(t — r') where Uq = 47th 2 a/M and a is the s-wave 
scattering length for elastic atom collisions. The Gross- 
Pitaevskii equation for the order parameter, then has 



Following Ref. we assume that the transverse di- 
mension of the cloud is sufficiently small and the cor- 
responding time scale sufficiently short that the trans- 
verse profile of the particle density can adjust to the 
equilibrium form appropriate for the local atomic den- 
sity. With this approximation, the problem becomes 
one-dimensional, and the solitary pulse can be described 
by a local velocity, v(z), and a local density of particles 
per unit length, cr(z) — J \^(x, y, z)\ 2 dxdy |l5| . The or- 
der parameter can then be written as a simple product 
7] with \&(r, t) = f(z,t)g(x,y,cr), where g is the equi- 
librium wavefunction for the transverse profile. If g is 
chosen to be normalized, J \g\ 2 dxdy — 1, the equations 
above imply that |/| 2 = a. 

In this problem there are three qualitatively distinct 
physical regimes and two corresponding transitions sep- 
arating them. The first of these transitions involves a 
change in the transverse profile of the cloud. For weak 
transverse confinement, \g\ 2 can be calculated using the 
Thomas-Fermi approximation and is simply a parabola. 
For transverse confinement of a strength sufficient to en- 
sure that the chemical potential is much smaller than 
hu)±, \g\ 2 becomes Gaussian. The second transition in- 
volves a change in the longitudinal profile of solitary 
waves. For weak transverse confinement, the atoms form 
a condensate. In the limit of strong interactions, strong 
transverse confinement and small linear densities, how- 
ever, bosonic atoms behave in a certain sense like non- 
interacting fermions. This transition involves no change 
in the (Gaussian) transverse profile since the chemical 
potential of the gas is much smaller than hu>± in both 
cases. 

We first consider the transition in the transverse di- 
rection. The critical value of era for this transition is 
determined by the condition that the kinetic energy as- 
sociated with transverse motion, ft 2 /MR 2 , is compara- 
ble to the typical interaction energy, uUq, where R is 
the transverse width of the cloud and n is the typi- 
cal (three-dimensional) density. Since n ~ N/(ZttR 2 ), 
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where N/Z = a is the density per unit length, one sees 
that the value of a at the crossover is <t C) i ~ 1/a. It is 
interesting to note that both energy scales vary as 1/R 2 
with the consequence that a c \ is independent of the os- 
cillator length <2_i_ = (fi/AfcjjJ 1 / 2 . 

Under typical experimental conditions (i.e., where a 
single trap rather than a series of tubes is used), a a » 1. 
The cloud is thus in the Thomas-Fermi regime. The situ- 
ation changes, however, with greater transverse confine- 
ment. To see this, we recall that the cloud expands along 
the long axis of the trap when it is squeezed transversely. 
Since Z increases, a — N/Z decreases. 

In the limit a a <§; 1, |<?| 2 has a Gaussian form, |g| 2 = 

e~( p / a ^ 2 /( 7ra i)- As shown in Ref. 0, / then satisfies 
the equation 

ihd t f = -(h 2 /2M)d 2 J + hu x (l + 2a\f\ 2 )f. (2) 

We see from this equation that / oc e -^±( 1 + 2acr o)t as 
\z\ — ► oo, where Co is the background linear density. 
Rewriting Eq. J5J using the variable w = fe iu± ( 1 + 2ttff o) i , 
we obtain 

ihd t w = ~(h 2 /2M)d 2 z w + huj ± 2a(\w\ 2 - cr a )w. (3) 

Equation pj includes a familiar (i.e., quadratic) nonlin- 
ear term and leads to a speed of sound, ci, which satisfies 
the equation Mc\ = 2huj±aoa. Since ctq = noira^, we see 
that Mc\ = n a U a /2 0. 

As shown in Refs. JllElj the density associated with 
a solitary wave of velocity u has the form 



a(z)/(TQ — 1 = —cos 2 #/cosh 2 (z cos#/£), 



(4) 



where 9 — sin (m/ci) and £ = 2^ (no)- Here £(n ) is 
the coherence length for n = ao/(na\), so that £ = 
a^/(2cr a) 1 ^ 2 . Further, the dispersion relation connect- 
ing the energy and momentum of a solitary wave, £(V), 
is given parametrically as £/£o = (4\/2~/3) cos 3 9, where 
So = huj±(aoa) 1 ' 2 aoax, and V/Vo — nu/\u\ — 29 — sin26*, 
where Vq = <j$Ti. 

We now consider the second, longitudinal transition. 
As noted above, the gas approaches the so-called Tonks- 
Girardeau limit, in which the bosons to some extent 
behave like non-interacting fermions, as the transverse 
confinement is increased. The motion of the atoms in 
this limit is effectively one-dimensional since transverse 
degrees of freedom are frozen out when huj± is much 
larger than the chemical potential. In the crossover re- 
gion, each atom occupies a length of order 1/cxo along 
the axis of the trap. The corresponding kinetic energy 
is on the order of Ti Oq/M. The typical interaction en- 
ergy, on the other hand, is on the order of nollo, where 
no ~ cr o j_ so that hqUo ~ h 2 aoa/ '(Maj_). Thus, the 
ratio between the interaction energy and the kinetic en- 
ergy is ~ a/ (aoo^) ~ [&o£,( n o)]~ 2 ■ Given the assumption 
of strong interactions, this quantity is much larger than 



unity for small values of Co (i- e -! for low densities), small 
values of a±_ (i.e., for strong transverse confinement), or 
large values of a (i.e., for strong interactions). The transi- 
tion between the two regimes takes place for a r ? ~ a/a 2 ^. 
For values of <7o much smaller than cr Cj 2 |lq|. 

ihdtf = (h 2 /2M)(-d 2 f + n 2 \f\*f). (5) 

In this case we see that / oc e -*0 2CT o ?i / 2M ) t as \ z \ — > 
oo. We thus rewrite Eq. © using the variable w — 

j- e i(rr 2 a 2 h/2M)t tQ obtain 

ifidtw = (h 2 /2M)[-d 2 w + tt 2 (H 4 - <r 2 )w}. (6) 

The speed of sound C2 is now given as Mc 2 , = n 2 h 2 ao/M 
or C2 = irhao/M. This is precisely the Fermi velocity in 
one dimension with kp = (Tqtt. 

As shown in Ref. Q , the density of the cloud associated 
with the solitary wave is now 



(t(z)/<j - 1 = -- 



3 cos 2 9 



2 + (1 + 3 sin 2 6») 1 /2 C osh(27rcr z cos 6) ' 

(7) 

where 9 = sin 1 (u/c2). The dispersion relation can again 
be expressed parametrically using 



E /E Q = ^l C0S 2 9\n 



2 + 3 cos 2 9 



(1 + 3 sin 2 0)V2 



with £ = % 2 <jI/M, and 



VjVo = -ir(£/£ ) tan^+cos -1 



3 sin 2 9 - 1 



(1 + 3 sin 2 6»)!/2 



(8) 



(9) 



It is natural to ask how one actually proceeds from 
one region to another. The transverse transition was 
studied in some detail in Refs. 0, 0- Here, we focus 
on the second, longitudinal transition using a formalism 
of more general applicability. The basic ingredient re- 
quired is the equation of state, i.e., the energy of the gas 
per unit length, e(a), as function of a. This can be calcu- 
lated from the Lieb-Liniger model [l7| and can be written 
as e(cr) = (h 2 a 3 /2M) §(7), where §(7) is a numerically 
known function and 7 = 2a/(<ra 2 L ). 

Given e(er), one can immediately determine the sound 
velocity as Ale 2 — a^d 2 € / da 2 (ao) . The energy of the 
solitary wave is 



£ = 



( h 2 dw* dw .. . \ , 

{2M^j; + ^-^^ ) + c ) dz (10) 



with e' — de/da. The final term in this equation rep- 
resents the energy of the background density of atoms 
and ensures convergence of the integral. Thus, C = 

-e(<r ) + eroe'(cTo). 

Equation (|lCj|) implies that w satisfies the equation 



ifidtw 



-(H 2 /2M)d 2 w + [e'(a) - e'(a )]w. (11) 
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Writing w — y/a e 1 ^ and separating the real and imagi- 
nary parts of Eq. (|llfl , we obtain the two hydrodynamic 
equations 

-Mu* {a - ao)2 (12) 

2(7 

and v = (h/M)d(j)/dz = u(l — cr /cr). Here, we have 
imposed the boundary condition that v — > as a — ► a§. 
Eq. I jlOjl can be written as 



£ = 



2 *2 



-(a-a )e'(a )]dz. (13) 
Combining the above equations, we obtain 

£ = 21 [e(<r) - e(a ) -{a- a )e'(a )] dz. (14) 



Finally the momentum of the solitary wave is given as 

(g = ^ dz. (15) 



V =M (a- <T )v(z) dz = Mu 



We thus see that, given knowledge of the energy per 
unit length e(cr), Eq. (|12|l allows us to determine the 
shape of the solitary wave <j{z) for a given velocity, it. 
Equations (|14f> and (|15l) then give £(u) and V(u), which 
can be combined to establish the dispersion relation, 
£ — £('P). While the profiles of the solitary wave and the 
corresponding dispersion relations are known in all three 
limits examined earlier H US ElEl E3. E q s - 03> " 
(|15|l allow us to interpolate between them. As indicated, 
our approach is quite general and merely requires knowl- 
edge of e(a). 

Let us now consider a specific example for the case 
7o = 2a/(<jQa\) = 1. It is convenient to express e 
explicitly in terms of a with the result that e(cr) = 
[Ti 2 Oq/2M] e(a/ao), where a reliable interpolation for- 
mula for e is given as e(y) = (7t 2 7q?/ 3 /3 + njoy 5 )/(lo + 
4 7o2/ + 47oy 2 + ny 3 ) with y = a(z)/a and re « 6.879. 

For 70 = 1, Eq. ((T^l) then has the form 



(dy/yV = ttV/3 + K y 5 
V dz J 1 + Ay + 4y 2 + ny 3 



(y-1) de/dy\ t 



(Mc 
\ha 



tt 2 /3 + re 
9 + re 

u 2 (y-1) 2 
c 2 y 



(16) 



where z = o$z. Given that de(y)/dy ~ 1.471 and 
d 2 e(y) / dy 2 » 1.872 for y = 1, we find that the sound 
velocity is given as (Mc/Tujq) 2 rs 0.936. For a given 
value of m/c, Eq. iflfijl gives the profile of the solitary wave 
<t(z)/<7o- Figure 1 shows such solutions for u/c = 0.1 
(bottom curve at z = 0), 0.4, 0.7, and 0.9 (top). 
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FIG. 1: The density <t(z)/<tq associated with a solitary wave 
given by the solution of Eq. 1161 . for 70 = 2a/(<roa^) = 1 and 
u/c = 0.1 (lowest graph at z = 0), 0.4, 0.7, and 0.9 (highest), 
where z = ooz. 



The energy and the momentum can then be calculated 
from Eqs. {Tl} and JTSJ: 



£(u) 



and 



M 



7r 2 y 3 /3 + rey° 



l+4y- 



4?/ 2 + rey 3 
-1.471 (y 



7r 2 /3 + re 

9 + re 
-l)]d~z, (17) 



= Mu / dz(y-iy/y. 



(18) 



The dispersion relation £ = £('P), initially obtained by 
Lieb, is shown as the continuous (lower) line in Fig. 2. 
The maximum value of V corresponds to a dark solitary 
wave with u = 0. For long wavelengths, the solitary 
waves become ordinary sound waves with the usual linear 
dispersion relation |6j, llOl llll |l2( . From Eq. I|ll|) one can 
also calculate the dispersion relation appropriate for the 
Bogoliubov mode. This has the usual form 



£ = y \r 2 /2M) 2 + {cvy 



(19) 



The Bogoliubov dispersion relation, £ = £(V), is shown 
as the dashed (higher) line in Fig. 2. 

The two conditions derived above establishing the typ- 
ical values of a that characterize the three regimes, 
fc.ia ~ 1 and (T Cj 2Ctj_ ~ a, give values of a which dif- 
fer by roughly one order of magnitude when the trans- 
verse confinement is strong as in Ref. For example, 
for a — 100 A and u±_ — 2ir x 70.7 kHz (which im- 
plies that a± ~ 0.04 /im), we see that er Cj i ~ 10 6 /cm and 
c c ,2 ~ 10 5 /cm. 

Another remarkable feature of this problem is that it 
is independent of the scattering length in the Tonks- 
Girardeau limit. Specifically, in the two limiting cases 
(i.e., the one-dimensional Gross-Pitaevskii equation and 
the Tonks- Girardeau equation), the characteristic widths 



4 





o 
o 




2.5 




2 






a> 


1.5 


d 




LU 


1 



0.5 ■ ^ ' 

^— ■ ■ ■ ■ ■ 

0.5 1 1.5 2 2.5 3 3.5 

Momentum 

FIG. 2: The dispersion relation £ = £(V) for solitary waves 
(lower) and for the Bogoliubov mode (higher). The momen- 
tum is measured in units of Tioo, and the energy is measured 
in units of h 2 ao/M. Here 70 = 2a/(aoa 2 L ) — 1. 
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FIG. 3: The profile of the solitary waves for u/c = 1/2 and 
three values of 70 = 2a/ (aoa^). In the lowest one (I) 70 = 
1/10, corresponding to the one-dimensional Gross-Pitaevskii 
limit, Eq. (J1J, in the middle one (II) 70 — 1, corresponding to 
the intermediate regime, Eq. I|16[l . and in the highest one (III) 
70 = 10, corresponding to the Tonks-Girardeau limit, Eq. 0. 



of solitary waves is on the order of the coherence length 
£ ~ aj^/((Toa) 1 / 2 and on the order of the interparticle 
spacing l/erg, respectively. This allows us to under- 
stand the qualitative features of the profiles shown in 
Fig. 3. For the largest value of ctq — 10er c , which ex- 
plores the Gross-Pitaevskii limit, the characteristic width 
of the disturbance (in units of (Jq 1 ) is on the order of 
(croa^/a) 1 / 2 ~ 1/^/to, which is larger than unity. In 
the opposite Tonks-Girardeau limit with a = ct c /10, the 
width is of order unity in units of a^ 1 . Therefore, as one 
moves from one regime to the other, the width of the soli- 
tary wave in the dimensionless units of Fig. 3 decreases 
and eventually reaches a constant value of order unity. 
This fact can be used as an experimental signature of 
the transition from one limit to the other. 

The present model is valid provided that the solitary 



waves have a width which is larger than the atom-atom 
spacing. While our predictions are thus expected to be 
reliable for waves whose size is larger than l/oo (e.g., 
sound waves), they are pushed to the limits of their va- 
lidity in the extreme cases of the Tonks-Girardeau limit 
(7 ^> 1) and for narrow (and thus slow) pulses. Predic- 
tions in these regions are expected to be qualitatively, 
but not quantitatively, correct Q . The Boboliubov spec- 
trum, on the other hand, agrees with the exact result 
for V — > and V — > 00 for all values of 7. 

Finally, as seen in Fig. 2, the dispersion relations for 
the two modes converge in the limit of long wavelength 
since both correspond to small amplitude sound waves 
in this limit. For shorter wavelengths, however, they di- 
verge. The characteristic momentum, V c , for which dif- 
ferences appear is on the order of h/£. Thus, Vc/'Po ~ 
■^/a/(aoa 2 L ) ~ i/70) which is on the order of unity at the 
transition and larger than unity as one approaches the 
Tonks-Girardeau regime. This behavior is quite different 
from that between the three-dimensional and the one- 
dimensional problems studied in Refs. 0,0,0, where 
the two modes have different energies and momenta. In 
this respect the present problem more closely resembles 
the homogeneous problem initially studied by Lieb 0, |(J . 
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